The entropy of hot QCD at large Nj: 
Successfuhy testing weak couphng techniques 

Jean-Paul Blaizot and Andreas Ipp 

ECT*, Villa Tambosi, Strada delle Tabarelle 286, 
1-38050 Villazzano Trento, Italy 

Anton Rebhan and Urko Reinosa 

Institut fiir Theoretische Physik, Technische Universitat Wien, 
Wiedner Hauptstr. 8-10, A-1040 Vienna, Austria 

Abstract 

It has been known for some time that the entropy of hot QCD is well 
reproduced by weak coupling techniques. These do not identify with per- 
turbation theory, known to have poor convergence properties, but involve 
resumming the physics of hard thermal loops in the non-perturbative 
derivable two-loop approximation. We test this approximation scheme in 
the limit of large flavor number (Nf) where the exact result can be calcu- 
lated for a wide range of couplings for which the influence of the Landau 
pole is negligible. Using full momentum-dependent fermionic self-energies 
instead of the previously only known weighted average value to order g^, 
the exact result for the entropy can be remarkably well reproduced by the 
HTL resummed theory for a natural choice of the renormalization scale, 
and this up to large values of the coupling. This gives confidence in the re- 
liability of weak coupling techniques when applied to the thermodynamics 
of QCD, at least for temperatures T > 3Tc. 

Introduction 

The entropy of hot QCD is well reproduced for temperatures T > 3Tc by the 
non-perturbative (f>-derivable two- loop approximation Its success has been 
established by comparison to lattice simulations Here we want to present 
an independent test jSj of this approximation scheme in the limit of large flavor 
number Nf 0]. 

Approximations based on the $ functional are constructed from the two- 
particle-irreducible (2PI) skeleton expansion where the thermodynamic po- 
tential is expressed in terms of dressed propagators according to 

n[G,S]/T = ^Tr logG-1 - iTiHG 

-Trlog5-i+ TrS5 + $[G,S'], (1) 

where "Tr" refers to full functional traces and $[G, S\ is the sum of 2PI "skele- 
ton" diagrams. The self-energies n[G, S\ G^^ -Gg ^ and I][G, S\ = S-^-S^^ 




Figure 1: Comparison of the averaged asymptotic thermal quark mass squared, 
M^, in various approximations. In this and the foUowing plots, the renormal- 
ization scale /2ms is varied around the FAC-m scale by factors of 2. 

are determined by the stationarity property S^l[G, S]/dG = = 5il.[G , S] / 6 S . 
For the entropy S = {dP/dT)^, with P ^ -Vl/V one obtains [J: 

S ^ f ^!!iMrinilogG-i- ImnRcGl 

Jk oT ^ 

-2tr / [In^ log^"' - ImEReS] . (2) 

In the large Nf limit for QCD the effective coupling g^^j = g'^Nf/2 
is kept finite as the limit Nf ^ cx) is taken. The thermodynamic potential 
to order iV° and to all orders in ^eff is then given by a set of ring-diagrams 
which can iDe resummed. Large-A^y QCD is no longer asymptotically free. Its 
renormalization scale dependence is determined (non-perturbatively) by the one- 
loop beta function pi^g^^) = g^jj/(67r^) which implies a Landau singularity. The 
two-loop entropy formula (O reduces in the large Nf limit to 

S~So = -tr / [inilogG-VGo'- IninReG] 

-2tr|^^^[ReEIm5o]. (3) 

Using the 1-loop self-energy 11, this formula actually reproduces the large Nf 
entropy to all orders in gcs_. The fermionic self-energy S from this formula only 
has to be evaluated on the light-cone. The bare self-energy Sb is given by 

EbW = Sth + Sb.vac = -5'CF^ l^'SoiQ + iqYG^.iQ). (4) 

The thermal piece Sth of the go integration is UV finite, while the real part 
of the vacuum piece ReEb.vac is in general logarithmically divergent, but finite 
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Figure 2: Asymptotic thermal quark mass squared as a function of k/T for 
g^ff{TTT) = 9. The exact large- A^/ result is compared to the NLO-HTL calcula- 
tion and its relevant average value. 

on the light-cone. Because of the Landau pole Al, we have to introduce a Eu- 
clidean invariant cutoff < A|^ , and the analytic continuation of Eb, vac to the 
light-cone requires a deformation of the integration path The numerically 
demanding calculations have been performed on the ECT* Teraflop cluster. 

HTL approximation of the self-consistent entropy 

The HTL effective action |H] is an effective action for soft modes with energy 
scales ^ gT. Inserting HTL into the entropy formula (O gives the correct 
contribution, but only a fraction of the plasmon term ~ g^. The larger part of 
the plasmon term arises from hard momentum scales where HTL is no longer 
accurate, namely from corrections to the leading-order asymptotic masses 

HT(fco = fc)-™L= 2^D, E±(fco = ±fc)-^ = — . (5) 

The NLO corrections to the asymptotic thermal masses are nontrivial functions 
of the momentum, 

Sml^{k) = Re(5HT(/so = k), = Re2kST,+ {ko = k), (6) 

and involve only a single HTL propagator and no HTL vertices for hard external 
momentum. These expressions can be evaluated only numerically (NLO-HTL). 
However, the weighted averages (NLA) are determined to order g^ 

As one can see in Fig. ^ the NLO-HTL-resummed result represents a con- 
siderable improvement over the HTL-resummed result truncated at order g^g 
for g^g > 4, which corresponds to rhn/T > 1. Figure [3 shows the results of a 
numerical calculation of the asymptotic thermal quark mass squared M^{k) = 
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Figure 3: Entropy in the large- iV/ limit, comparing the exact large- A'^y result to 
the strictly perturbative expansion through order g\,^^ the NLA result, and the 
NLO-HTL result. 

2A:ReI](fco ^ k) for ffeff (mms = ttT) = 9, normalized by T'^Cf/Nf. The exact 
(nonperturbative) result obtained in the large- A^y limit is given by the full lines. 
In Fig. O the exact and the NLO-HTL results for the entropy are compared to 
simpler approximations. It is remarkable that a full HTL resummation using 
momentum-dependent asymptotic quark masses can extend the validity of the 
approximation in the large Nf limit up to rhu/T ~ 2.5. For the entropy of (small 
Nf) QCD we therefore expect likewise improvement once full HTL-resummed 
momentum-dependent asymptotic gluon masses have been implemented. 
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